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Abstract 

Ba2CoGe207, crystallizing in the noncentrosymmetric but nonpolar P 42 im structure, belongs 
to a special class of multiferroic materials, whose properties are predetermined by the availability of 
the rotoinversion symmetry. Unlike inversion, the rotoinversion symmetry can be easily destroyed 
by the magnetization. Moreover, due to specific structural pattern in the xy plane, in which the 
magnetic Co^''" ions are separated by the nonmagnetic Ge04 tetraherda, the magnetic structure of 
Ba2CoGe207 is relatively soft. Altogether, this leads to the rich variety of multiferroic properties of 
Ba2CoGe207, where the magnetic structure can be easily deformed by the magnetic field, inducing 
the net electric polarization in the direction depending on the magnetic symmetry of the system, 
which itself depends on the direction of the magnetic field. In this paper, we show that all these 
properties can be successfully explained on the basis of realistic low-energy model, derived from 
the first-principles electronic structure calculations for the magnetically active Co 3 d bands, and 
the Berry-phase theory of electronic polarization. Particularly, we argue that the magnetization 
induced electric polarization in Ba2CoGe207 is essentially local and expressed via the expectation 
values (p) = Tr[p'D] of some dipole matrices p, calculated in the Wannier basis of the model, 
and the site-diagonal density matrices V of the magnetic Co sites. Thus, the basic aspects of the 
behavior of Ba2CoGe207 can be understood already in the atomic limit, where both magnetic 
anisotropy and magnetoelectric coupling are specified by V. Then, the macroscopic polarization 
can be found as a superposition of (p) of the individual Co sites. We discuss the behavior of 
interatomic magnetic interactions, main contributions to the magnetocrystalline anisotropy and 
the spin canting in the xy plane, as well as the similarities and differences of the proposed picture 
from the phenomenological model of spin-dependent p-d hybridization. 
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I. INTRODUCTION 


Magnetically driven ferroelectricity is one of the major topics in the condensed matter 
physics today. Even though the basic crystallographic symmetry of a magnetic material 
may not allow for the existence of a spontaneous electric polarization, this symmetry can be 
further lowered by the magnetic order, which in some cases gives rise to the ferroelectric (FE) 
activity. After discovery of such effect in TbMnOs,- where the inversion symmetry is broken 
by some complex noncollinear magnetic order, there is a large number of experimental and 
theoretical proposals pointing at the existence of a similar effect in other magnetic materials, 
which are called multiferroics.- 


Ba2CoGe207 and related compounds take a special place among multiferroic materials. 
They crystallize in the noncentrosymmetric P42im structure,- which nonetheless does not 
permit the FE effect because of the fourfold rotoinversion symmetry. Nevertheless, unlike 
inversion, the rotoinversion symmetry can be relatively easy destroyed by the magnetic order, 
by canting the magnetic moments out of the rotoinversion axis. Thus, one unique aspect of 
Ba2CoGe207 is that the ferroelectricity in this compound can be induced relatively simple C- 
type antiferromagnetic (AFM) order. Furthermore, the magnetic structure of Ba2CoGe207 
is relatively soft: It consists of the C0O4 tetrahedra, which are interconnected with each 
other via nonmagnetic Ge04 tetrahedra (see Fig. \T^). Since two magnetic Co sites in this 
structure are separated by the long Co-O-Ge-O-Co paths, the exchange interactions between 
them are expected to be small and the magnetic structure can be easily deformed by the 
external magnetic field. The possibilities of easy manipulation by the magnetic structure 
and switching the FE properties have attracted a great deal of attention to Ba2CoGe207, 
and today it was demonstrated in many details how the electric polarization in Ba2CoGe207 
can be tuned by the magnetic field (Refs. 4 andO) as well as the uniaxial stress (Ref. 1^. 


The experimental behavior of the electric polarization in Ba2CoGe207 is frequently inter¬ 
preted basing on the model model of spin-dependent p-d hybridization: if e is the direction 
of the spin magnetic moment at the Co site, located in the origin, and v is the position 
of a ligand oxygen atom, relative to this origin, the bond Co-0 contributes to the local 
electric dipole moment as p oc (e ■ This simple phenomenological expression is able 

to capture the symmetry properties of the electric polarization and, thus, explain the be¬ 
havior of this polarization on a phenomenological level. It is therefore not surprising that a 
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FIG. 1 . (Color online) Crystal and spin magnetic structure of Ba2CoGe207: (a) General view 
on the crystal structure in the xy plane. The Co and O atoms are indicated by the medium red 
and yellow spheres, respectively, the Ba atoms are indicated by the big green spheres, and the 
Ge atoms are indicated by the small blue spheres. The Mn04 and Ge04 tetrahedra are colored 
red and green, respectively, (b) Details of the rotations of Mn04 tetrahedra associated with two 
Mn sites in the xy plane, (c) and (d) Two possible spin magnetic structures, realized without 
external magnetic field. The light (cyan) arrows indicate the directions of spins favored by the 
single-ion anisotropy, while the dark (blue) arrows are the joint effect of the single-ion anisotropy 
and interatomic exchange interactions. 

qualitatively similar behavior of the electric polarization was obtained in the first-principles 
electronic structure calculations,- which generally support the model of spin-dependent p-d 
hybridization. 

On the other hand, the correct quantum-mechanical definition of the electric polarization 
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is solids should be based on the Berry-phase theory.— This theory provides not only an effi¬ 
cient computational framework, which is used today in most of the first-principles electronic 
structure calculations, but also appears to be a good starting point for the construction 
of realistic microscopic models, explaining the origin and basic aspects of the behavior of 
electric polarization in various types of compounds.—-— In this work we will pursue this 
strategy for the analysis of electric polarization in Ba2CoGe207. First, we will show that, 
in the atomic (Wannier) basis, the electric polarization consists of two parts: the local one, 
which is expressed via the expectation value (p) = TrlpI)] of some local dipole matrices 
p and site-diagonal density matrices P, and the anomalous one, which is solely related to 
the phases of the coefficients of the expansion of the wavefunctions in the basis of atomic 
Wannier orbitals. Similar classification holds for the orbital magnetization.— Then, if the 
crystal structure possesses the inversion symmetry, the local term becomes inactive and the 
electric polarization, induced by the magnetic inversion symmetry breaking, is anomalous in 
origin. Such situation occurs, for instance, in multiferroic manganites with the orthorhombic 
structure.— In Ba2CoGe207, however, the situation is exactly the opposite: to a good ap¬ 
proximation, the electric polarization is expressed as a sum of local electric dipole moments 
of individual Co sites, while the anomalous contribution is negligibly small. This explains 
many experimental details of the behavior of electric polarization in Ba2CoGe207 as well as 
basic difference of this compound from other multiferroic materials. 

The rest of the paper is organized as follows. The method is described in Sec. [TTl We 
use the same strategy as in our previous works, devoted to the analysis of multiferroic 
properties of various transition-metal oxides on the basis of effective realistic low-energy 
models, derived from the first-principles calculations.— Therefore, in Sec. Ill Al we will 
remind the main details of the construction of such effective low-energy model. In Sec. Ill Bl 
we will analyze the main contributions to the electric polarization in the case of basis, when 
the Bloch wavefunction is expanded over some atomic Wannier orbitals, centered at the 
magnetic sites. Since we do not consider the magnetostriction effects, which can low the 
original P 42 im symmetry of the crystal structure, the considered electric polarization is in 
fact the electronic one.- In Sec. lIIIl we will present results of our calculations for Ba2CoGe207, 
starting from the atomic limit and, then, consecutively considering the effect of interatomic 
exchange interactions and the external magnetic field. Finally, in Sec. lIVi we will summarize 
our work. In the Appendix we will estimate the ferromagnetic (FM) contribution to the 
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interatomic exchange interaction cause by the magnetic polarization of the oxygen 2 p band 
and discuss how such effect can be evaluated starting from the low-energy electron model 
for the magnetic Co 3 d bands. 


II. METHOD 

A. Effective low-energy model 

In this section, we briefly remind the reader the main ideas behind the construction of the 
effective low-energy electron model. The details can be found in the review article fRef. 1 151 ). 
The model Hamiltonian, 

ij a 0 i a/S'yS 

is formulated in the basis of Wannier orbitals { 4 >ia}, which are constructed for the mag¬ 
netically active Co 3 d bands near the Fermi level, starting from the band structure in the 
local-density approximation (LDA) without spin-orbit (SO) coupling (Fig. ( 21 ). Here, each 



Energy (eV) 

FIG. 2 . (Color online) Total and partial densities of states of Ba2CoGe207 in the local density 
approximation. The shaded light (blue) area shows the contribution of the Co 3 d states. The 
positions of the main bands are indicated by symbols. The Fermi level is at zero energy (shown by 
dot-dashed line). 
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Greek symbol (a, j5^ 7, or J) stands for the combination of spin (cr^, ^7/3, cr^, or gs) and 
orbital (a, 6, c, or d) indices, for which we adapt the following order: xy^ yz, 3 ^^— zx, 
x'^—y‘^. Each lattice point i (j) is specified by the position r (r') of the atomic site in the 
primitive cell and the lattice translation R. Hence, the basis orbitals ^iQ;(r) = R— r) 

are centered in the lattice point (R+r) and labeled by the indices r and a. Moreover, they 
satisfies the orthonormality condition: 

{4>T'a'{r - R' - T')|0™(r - R - t)) = (2) 


In our case, the Wannier function were calculated using the projector-operator method 


(Refs 




and 


161 ) and orthonormal linear muffin-tin orbitals (LMTO’s, Ref. 


171 ) as the 


trial wave functions. Typically such procedure allows us to generate well localized Wannier 
functions, that is guaranteed by the good localization of LMTO’s themselves. Then, the 
one-electron part of the model ([T]) is identified with the matrix elements of LDA Hamilto¬ 
nian ("Hloa) in the Wannier basis: = (0TQ(i'“'^)|’HLDA|^/^T'y3(i'—R— t^)). Since the 

Wannier basis is complete in the low-energy part of the spectrum, the construction is exact 
in the sense that the band structure, obtained from exactly coincides with the one 

of LDA. 

All Wannier basis and one-electron parameters were first computed without the 

SO interaction. Then, the SO interaction at each atomic site was included in the “second- 
variation step”, in the basis of the nonrelativistic Wannier functions: A7{so|0i/3(i’))i 

as explained in Ref. [ 15 . 

Matrix elements of screened Coulomb interactions can be also computed in the Wannier 
basis as 

Uap-yS = j dr j drV*^(r)(/)i^(r)'yscr(r,r')0;y(r')(^w(r'), 

where Vscri'J^, r') is obtained using the constrained RPA technique.— Then, Uscr(r, r') does not 
depend on spin variables and Since RPA is very time consuming 

technique, we employ additional simplifications, which were discussed in Ref. Il 5 |. Namely, 
first we evaluate the screened Coulomb and exchange interactions between atomic Co 3d 
orbitals, using fast and more suitable for these purposes constrained LDA technique. After 
that, we consider additional channel of screening caused by the 3d —>■ 3d transitions in 
the polarization function of constrained RPA and projecting this function onto atomic 3 d 
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orbitals. The so obtained parameters of screened Coulomb interactions are well consistent 
with results of full-scale constrained RPA calculations without additional simplifications. 

All calculations have been performed for the room temperature PA 2 im structure re¬ 
ported in Ref. The corresponding parameters of the low-energy model are summarized in 
supplemental materials.— 

After the construction, the model is solved in the unrestricted Hartree-Fock (HF) approx¬ 
imation, which is well justified for the considered case where the degeneracy of the ground 
state is lifted by the crystal distortion.— 


B. Electronic polarization 


The electronic polarization can be computed in the reciprocal space, using the formula 
of King-Smith and Vanderbilt:- 


P = 


/o ^ ^ (^nk| 5 

[27:Y ^ Jbz 


( 3 ) 


where Wnk(i*) = is the cell-periodic eigenstate of the Hamiltonian = 

g-ikr^gikr^ the Summation runs over the occupied bands (n), the k-space integration goes 
over the first Brillouin zone (BZ), and —e (e > 0 ) is the electron charge. In our case, each 
' 0 nk(i') is expanded in the basis of Wannier orbitals 0 T-c^(r —R —r), used for the construction 
of the low-energy model: 


1 


^„k(r) = ^ E - R - 


( 4 ) 


R,tq 


where N is the number of primitive cells. Then, the k-space gradient of w^k will have two 
contributions: 


SkWnk = -^^(r-R - - R - r) 

^ ^ Rtq 

+ E - R ~ t), 


( 5 ) 


Rra 


and the electronic polarization P will also include two terms: 


dk 


P — Q (CnklPkl'^nk) ^ (^nkl<5kCnk) ■ 


' BZ 


le 

(27r)^ 


( 6 ) 


' BZ 





Here, |cnk) denotes the column vector |cnk) = [c^k]- Pk matrix 

where 

Pk“'™ Y.(^r'a'{r - r')|r|,^™(r - R - 


( 7 ) 


R 


V is the unit cell volume, and Q = ( 27 r)^/V is the first BZ volume. Moreover, since P is 
understood as the change of the polarization in the process of adiabatic symmetry lowering 
(Ref. 0 ), which in our case is driven by the magnetic degrees of freedom, here and below we 
drop all nonmagnetic contributions to ^ which are irrelevant to the magnetic symmetry 
lowering. 

If the transition-metal site is located in the inversion center, the first term of Eq. © iden¬ 
tically vanishes. Such situation is realized, for instance, in multiferroic manganites, crystal¬ 
lizing in the orthorhombic Pbnm structure, where the FE activity is entirely “anomalous” 
and associated with the second term of Eq. ©.— However, in Ba2CoGe207 the situation 
appears to be exactly the opposite: the first term dominates, while the second contribution 
is negligibly small (about 0.1 pC/m^ in the magnetic ground state). Furthermore, due to 
the orthogonality condition (| 2 ]), the leading contribution to pj^“ comes from the Wannier 
functions centered at the same atomic site. Then, one can impose in Eq. ([Tj) the additional 
condition t' = R+r, which yields ~ p^'^Sr^r^ where p"^^ does not depend on k. 

Thus, the electronic polarization in Ba2CoGe207 will be given by the sum 


of the local electric dipoles: 




Pr = (Pr) = 


( 8 ) 


( 9 ) 


where Tr^ is the trace over a. Each such dipole is given by the expectation value (p^) of the 
dipole matrix p^ = [Pt^] and the spin-independent part of the density matrix T)^ = 
at the site r: 


-'nk’-'nk 


In order to evaluate p^ we use the LMTO method and expand each (/>ra(i’) is the basis 
{Xv^} of linear muffin-tin orbitals, which can be viewed as orthonormalized atomic-like 
Wannier orbitals, constructed in the whole region of valence states.— Moreover, we shift 
each site r to the origin, that does not affect the polarization change caused by the magnetic 
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symmetry lowering. Then, without SO interaction, we have 


'^ra(r) = IraXvbir -V + t). (10) 

vb 

where, v = r corresponds to the “head” of the Wannier function, centered at the Co-site 
r, while all other contributions describe the “tails” of the Wannier functions, spreading to 
the neighboring sites v. Then, around each such site we identically present the position 
operator as r = (v—t) + (r— i;+t) and assume that the leading contribution to comes 
from the first term. This is reasonable because, at each site, (r — + r) couples the atomic 

states with different parity, which are typically well separated in energy. Then, the matrix 
elements can be easily evaluated as 


= -V 


vb 


Then, by considering the leading contributions of four oxygen sites surrounding each Co 
atom, one obtains the following matrices (in /iC/m^): 


Pl,2x = 


Pl,2y - 


and 


Pl,2z = 



0 

t3509 

0 

6357 

0 


t3509 

0 

-801 

0 

-3553 


0 

-801 

0 

±303 

0 



6357 

0 

±303 

0 

=f6142 

V 

0 

-3553 

0 

=f6142 

0 


/ 

0 

6357 

0 

±3509 

0 

1 


6357 

0 

=f303 

0 

t6142 


0 

=f303 

0 

-801 

0 


±3509 

0 

-801 

0 

3553 

V 

0 

=f6142 

0 

3553 

0 

) 


0 

0 

-3682 

0 

0 

\ 


0 

±3394 

0 

3524 

0 


- 

3682 

0 

0 

0 

±2870 


0 

3524 

0 

T3394 

0 


v 

0 

0 

±2870 

0 

0 

/ 


( 11 ) 


( 12 ) 


(13) 
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where the upper (lower) signs correspond to the Co-sites 1 ( 2 ) in Fig. [H As expected, the 
matrices — {Pi,2x^Pi,2y^Pi,2z) obey the P 42 im symmetry. If the density matrix Pi^2 
obeys the same symmetry, all local electric dipoles will vanish and there will be no net 
polarization in the ground state. However, if the symmetry of Vi^2 is lowered by magnetism, 
one can expect the appearance of local electric dipoles with some order between sites 1 and 
2 . If this order is antiferroelectric, there will be no net polarization. Nevertheless, if this 
order permits a FE component, the system will exhibit a finite net polarization. The details 
of such symmetry lowering will depend on directions of magnetic moments at the sites 1 and 
2 . For instance, if the spins are parallel to the 2: axis, the fourfold rotoinversion around 2:, 
5 |, will remain among the symmetry operations of the magnetic space group, and there will 
be no local electric dipoles. However, if the moments lie in the xy plane, the rotoinversion 
51 is replaced by the regular twofold rotation, ( 7 |, which allows for the existence of local 
electric dipoles parallel to 2:. For an arbitrary direction of spin, the symmetry will be further 
lowered, and the local electric dipoles may have all three components. However, the relative 
directions of dipoles at the sites 1 and 2 will depend on other symmetry operations. Below, 
we will investigate this magnetic symmetry lowering more in details. 

To conclude this section, we would like to emphasize that Eqs. ([H|) and Q are the cor¬ 
rect quantum-mechanical expressions for the electric polarization in Ba2CoGe2 07, which 
are based on the very general Berry-phase theory.- These expressions allows us clarify the 
microscopic origin of the polarization, which is frequently ascribed to the spin-dependent p-d 
hybridization.- The p-d hybridization does play a very important role in this material as it 
defines the matrix elements of in Eq. (Q). Moreover, the contribution of each Co -0 bond 
to the local electric dipole moment is indeed proportional to the vector (v—t), connecting 
the Co and O sites. However, the matrix elements themselves do not depend on the 
spin state. In this sense, the spin dependence of the hybridization is not the most impor¬ 
tant aspect. The spin dependence of the electric polarization in our picture comes from the 
matrix elements of the density matrix, which describes the effect of the SO coupling at the 
transition-metal sites. The local magnetic moment deforms the electron density around the 
Co sites. This deformation, which itself depends on the direction of local magnetic moment, 
spreads to the neighboring oxygen sites via the tails of the Wannier functions and produces 
finite electric moment. This is the basic microscopic picture underlying the formation of the 
local electric moments in Ba2CoGe207. 
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III. RESULTS AND DISCUSSIONS 


The first important question we need to address is the local properties developed at each 
of the Co sites. For these purposes we set all transfer integrals in Eq. ([T|) equal to zero 
and solve the model in the mean-field approximation separately for each Co site. The SO 
interaction in these calculations is treated in the frameworks of the self-consistent linear 
response (SCLR) theory.— For each direction of spin, it gives us the self-consistent HF 
potential in the first order of the SO coupling, which can be used again as the input of HF 
equations in order to obtain at the output the change of the density matrix and the total 
energy beyond the first order. We would like to emphasize that the considered below effects 
are beyond the first order of the SO coupling. 

The obtained polarization is in the good agreement with the regular self-consistent HF 
calculations, which can be performed for the high-symmetry points. For the in-plane ro¬ 
tations of the spin magnetization, the results are summarized in Fig. Q. In this case, 



FIG. 3 . (Color online) Electronic polarization, absolute value of orbital magnetization, and total 
energy depending on the direction e = (cos (^, sin 0 , 0 ) of spin magnetization in the xy plane. The 
lines show results of self-consistent linear response theory for the spin-orbit coupling/^S The symbols 
show results of unrestricted Hartree-Fock calculations for the high-symmetry points. 

and Py identically vanish, while Pz{(l>) exhibits the characteristic cosine-like behavior.- As 
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expected, each Pz{<i)) takes its minimum value at the angle 0^, which specifies the direction 
of the upper 0-0 bond of the C0O4 tetrahedron (see Fig. [H). Note also the minus sign 
in Eq. ([TD, which means that the considered polarization is purely electronic. Since two 
Co sites in Ba2CoGe207 are connected by the symmetry operation {CJ|cti/2+a2/2} (the 
twofold rotation around x, ( 7 f, associated with the half of the primitive translations in the 
xy plane), the angles 4>2 and (j)i satisfy the property: (j)^ = 7 rl 2 — (l)^, where the experimental 
value of (pi is 67 °.- The maximums of Pz{(p) are at cp = +7r/2, which specify the directions 

of the lower 0-0 bonds of the C0O4 tetrahedra. Correspondingly, Pz{(p) = 0 when the spin 
moment is aligned in between the upper and lower 0-0 bonds {(p = -h 7r/4), and the 
negative and positive contributions to the polarization cancel each other. 

Very generally, the single-ion (SI) anisotropy energy in tetragonal systems has the fol¬ 
lowing form:— 

Egi = Ki sin^ 6 K 2 sin^ 6 K 2 sin'^ 6 cos 4 ((^ — pP). 


By fitting the total energies, obtained in the SCLR calculations, one can find that Ki = 


— 2.307 meV, K2 = 0.072 meV, and K2 = — 0.018 meV. The obtained value o 
sistent with the experimental estimate Ki — 2.327 meV, reported in Refs. 


Ki is well con- 
22 and 23. The 


parameter K2, controlling the in-plane anisotropy, is small but finite. Anyway, this behavior 
is different from the conventional 5 = 3/2 spin-only model, where no in-plane anisotropy is 
expected.— The difference is caused by an appreciable orbital magnetization (~ 0.434 per 
Co site), which is unquenched in the xy plane. Moreover, this orbital magnetization exhibits 
a sizable anisotropy (~ 0.005 y-Q per Co site - see Fig. [ 3 |), depending on the direction in the 
xy plane. Then, it is reasonable to expect that the magnetocrystalline anisotropy energy 
should be related to the behavior of the orbital magnetization.— Indeed, we observe a close 
correlation between behavior of the SI anisotropy energy and the anisotropy of the orbital 
magnetization in the xy plane (see Fig. [ 31 ). Thus, similar to the polarization, the origin of 
the in-plane SI anisotropy is related to the orbital degrees of freedom, which become active 
due to the magnetic symmetry lowering. The SI anisotropy can be further controlled by 
applying the uniaxial stress .- 

By summarizing the behavior of the SI anisotropy energy, the spin moments are expected 
to lie in the xy plane and, at each Co site, be parallel to either upper or lower 0-0 bond. 
Therefore, as far as the SI anisotropy energy is concerned, one expects two magnetic config¬ 
urations, which are depicted in Figs. [Tt and [Hi. In the first case (Fig. [T]:), the spins at the 
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site 1 and 2 are parallel to, respectively, lower and upper 0-0 bonds and hence this phase 
if antiferroelectric. In the second case (Fig. [Ill), all spins are parallel to the upper bonds, 
giving rise to the FE order. In the atomic limit, these two configurations are degenerate. 

Another factor, controlling the relative direction of spins, is the interatomic exchange 
interactions. The isotropic part of these interactions can be computed by considering the 
infinitesimal rotations of spins, which provides the local mapping of the total energies onto 
the spin (Heisenberg) model ^ij^i ' with denoting the direction of spin 

at the site z.— There are six types of nonvanishing exchange interactions, which are ex¬ 
plained in Fig. m The values of these interactions, obtained in the C-type AFM state, are 



FIG. 4. (Color online) Main parameters of isotropic exchange interactions. 


Ji = 0.012 meV, J 2 = -2.730 meV, J 3 = -0.046 meV, J 4 = -0.071 meV, J 5 = -0.007 
meV, and Jq = —0.012 meV. Thus, the exchange interaction J 2 = J stabilizes the AFM 
coupling in the xy plane, while the FM coupling between the planes is stabilized by the com¬ 
bination of Ji, J3, and J5. The leading interaction J is relatively weak (compared to other 
transition-metal oxides), that is directly related with the crystal structure of Ba 2 CoGe 207 , 
where neighboring Co sites in the xy plane are separated by relatively long Co-O-Ge-O-Co 
paths. Nevertheless, the value of J seems to be overestimated in comparison with the exper¬ 
imental data. For example, the Neel temperature Tn ~ 34 K, evaluated using Tyablikov’s 
random-phase approximation,— is about five time larger than the experimental 7 n = 6.7 K.- 
Similar disagreement is found for the exchange coup ing itself: our calculations overestimate 


the experimental J, reported in Refs. [22, l23|, andl27|, by the same order of magnitude. This 
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seems to be a negative aspect of the low-energy model ([T]), which was constructed only for 
the Co 3d bands and neglects several important contributions to the magnetic properties of 
Ba 2 CoGe 207 , related to the magnetic polarization of the O 2p band. This is also the main 
reason why we are not able to obtain a good quantitative agreement with the experimental 
data for the behavior of electric polarization in the magnetic field: although our low-energy 
model correctly reproduces the main tendencies, the magnetic field should by additionally 
scaled (by the same factor as the exchange coupling J) to be compared with the experi¬ 
mental data. In the Appendix, we will evaluate the change of the magnetic energy, caused 
by the polarization of the O 2p band, and show that this mechanism, which favors the FM 
alignment, can indeed reduce the effective AFM coupling J. Nevertheless, this mechanism 
overestimates the FM contribution to J and alone does not resolve the quantitative dis¬ 
agreement with the experimental data, which should involve additional factors, such as the 
exchange striction below 

The Dzyaloshinskii-Moriya (DM) interactions can be calculated by applying SCLR theory 
for the SO coupling and considering a mixed perturbation, where the SO interaction is 
combined with rotations of the spin magnetization.— The DM interactions between nearest 
neighbors along 2 : are forbidden by the symmetry. Therefore, the strongest interaction is 
again ^2 = (d^,, dy, dz) , which takes place between nearest neighbors in the xy plane. For the 
bond, connecting the sites 1 and 2 in Fig. [Qi, this vector is given by d 2 = (—5, 5, —6) fieV. 
Similar parameters for other bonds can be obtained by applying the symmetry operations 
of the space group P42im, which will change the signs of d^ and dy. Thus, 6,2 is comparable 
with K 2 and can also contribute to the canting of spins. 

The AFM interaction J enforces the collinear alignment between spins in the xy plane 
(see Fig. d]), while the SI anisotropy and DM interactions result in a small canting of spins, 
which is given by:— 


5 ( 1 ) 


Ji 2 cos 4$ sin 405 + dz 
27 


(14) 


where $ = |(0i + ^ 2 ) is the average azimuthal angle, formed by the spins 1 and 2. It is 
interesting that the magnitude of the canting depends on <I>, which contributes to the SI 
anisotropy, but not to the DM energy. Indeed, for the antiferromagnetically coupled spins, 
being parallel to the [100] and [110] axes in the xy plane, cos44> is equal to 1 and —1, 
respectively. Therefore, in the first case (Fig. [Tt), the effects of the SI anisotropy and DM 
interactions will partly compensate each other (note that sin 405 < 0), while in the second 
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case (Fig. [Hi), these two terms will collaborate, that leads to larger spin canting. This 
analysis is totally consistent with results of unrestricted HF calculations for the low-energy 
model dU. 

The exchange coupling J, in the combination with the SI anisotropy 5 
generacy between states e||[100] and e||[110]. However, the corresponding energy difference 
(per one formula unit),— 

AF; = 27^^cos40?, 

is very small (about 1 fieV), mainly because cos4(/)5 is small. This is the main reason why 
the direction of spins in the xy plane cannot be not easily determined experimentally.— 

As was discussed in Sec. EB the finite value of the polarization is due to the magnetic 
symmetry lowering, which is reflected in the behavior of the density matrices Pi, 2 - The 
latter can be identically presented as P 12 = P? 2 + ^Pi ,25 where PJ 2 average density, 

obeying the P42im symmetry, and ( 5 Pi ,2 is a perturbation, which depends on the direction 
of spin e. Straightforward unrestricted HF calculations yield (in 10“^) 


and 


Re(^Pi,2 


/ 0 0 t3.92 0 

0 3.43 0 t1-90 

T3.92 0 0 0 

0 t1.90 0 -3.43 

\ 0 0 -1.30 0 


0 ^ 
0 

-1.30 

0 

0 / 


Re(^Pi,2 


/ 0 0 

0 ±1.85 

-0.16 0 
0 - 2.88 
\ 0 0 


-0.16 0 0 \ 

0 - 2.88 0 

0 0 t3.90 

0 t1-85 0 

T3.90 0 0 / 


(15) 


(16) 


for the solutions with the spin magnetization being parallel to the axes [100] and [110], 
respectively (see Figs. [TJd and c), where the upper (lower) signs correspond to the Co-sites 
1 (2) in Fig. [H By combining these matrices with pi, 2 z? given by Eq. (1131) . one obtains the 
following contributions of the sites 1 and 2 to the electric polarization: Pi, 2 z = ±31.2 yC/iv? 
and Pi ,22 = —28.9 /iC/m^ for e||[100] and e||[110], respectively. Thus, the magnetic structure 
with e||[100] preserves the symmetry operation {Cf |ai/ 2 ±a 2 / 2 }, connecting the sites 1 and 
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2, which results in the antiferroelectric state. On the other hand, the magnetic structure 
with e||[110] breaks this symmetry, giving rise to the FE order. The total polarization 
Pz = -P 12 +-P 22 = 57.8 fiC/iTp is in fair agreement with the experimental data. According to 
Eqs. 1131), (1T3I) . and such a behavior is related to the phases of the matrix elements of 
Re(5Pi^2? which interplay with the phases of pi, 22 ‘ The phases are organized in such a way 
that for e||[100] and e||[110] their interplay yields P 22 = Piz and P 22 = “R 12 , respectively. 
Note also that the obtained matrices Re5Pi^2 do not couple with pi^ 2 x and Pi, 2 y, so that the 
X and y components of the polarization are identically equal to zero. 

Next, we discuss how the electric polarization can be controlled by the external magnetic 
field H. First, we consider the situation, where the AFM spins are parallel to the [110] axis 
(Fig. (Hi) and apply H along the perpendicular direction [110]. The results are summarized 
in Fig. El where for an easier comparison with experimental data we plot — P 2 .— In this case. 



FIG. 5. (Color online) Electric polarization and net spin magnetization in the external magnetic 
field parallel to the [IlO] axis, as obtained from the solution of the low-energy electron model in 
the Hartree-Fock approximation. 

H controls the magnitude of the spin canting and the direction of spins relative to the upper 
(lower) 0-0 bonds. When both spins are parallel to the upper 0-0 bonds, —Pz takes the 
maximal value. The corresponding magnetic field can be easily found from the analysis of 
the spin Hamiltonian, which yields 

™ (^1 “ 0 • 
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Using above values of J and 0?, ^ ~ 2 and S = 3/2, Hm can be evaluated as 47 T, which is 
in the very good agreement with the results of HF calculations for the low-energy electron 
model displayed in Fig. 0 Nevertheless, is overestimated by factor five in comparison 
with the experimental data,-*^ following similar overestimation of J, as was explained above. 

When the spins are aligned in between the upper and lower 0-0 bonds, is equal to 
zero. The maximal value of the polarization at P = is Pm ~ 90 /iC/m^, which is in fair 
agreement with the experimental data (about 120 /iC/m^, Ref. 0). 

Since the in-plane anisotropy is small, the spins in the xy plane can be easily rotated by 
the magnetic field, which couples to the net magnetization. The results of such calculations 
are displayed in Fig. El In this case, the value of the magnetic field also plays an important 



FIG. 6. (Color online) Behavior of electric polarization (total and partial contributions of the 
Co sites 1 and 2) under the rotation of the external magnetic field in the plane xy, as obtained 
from the solution of the low-energy electron model in the Hartree-Fock approximation. Here, the 
angle 0 specifies the direction of the magnetic field, while the antiferromagnetic component of the 
magnetization is perpendicular to the filed. 

role, as it controls the angle between spins at two Co sublattices: if H ^ Hm, the individual 
contributions Pu and P22 change “in phase”, and for (j) = 7r/4 (modulo 7r/2) IP12+P22I 
reaches its maximal possible value Pm- However, if P ^ Pm, there is some “dephasing” and 
IP 12 +P 22 I is smaller than Pm- 

As was already discussed in Sec. Ill Bi for an arbitrary direction of spins, the original 
P42im symmetry is completely destroyed and the polarization P can also have an arbitrary 
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direction. Similar to the in-plane rotations (Fig. [3|), this behavior can be well understood 
already in the atomic limit, by considering the change of the density matrix, induced by 
the SO interaction at a given Co site, which couples with the electric dipole matrix p^. For 
these purposes we again employ the SCLR method. It is convenient to start with the AFM 
configuration of spins parallel to the [110] axis and rotate them out of the xy plane. The 
results of such calculations are summarized in Fig. [71 One can clearly see that, in addition 



FIG. 7. (Color online) Electronic polarization and total energy depending on the direction e = 
(cos 0 sin sin 0 sin cos0) of the spin magnetization, rotated out of the xy plane for 0 = 7 r /4 and 
(j) = 57 r/ 4 , as obtained in the self-consistent linear response theory for the spin-orbit coupling^^ 
The azimuthal angles 0 = 7 r /4 and (p = 57 r /4 at the sites 1 and 2, respectively, were chosen to 
simulate the antiferromagnetic spin alignment in the xy plane. 

to Pz, there are finite perpendicular components of the polarization, and which obey 
certain symmetry rules and replicate the geometry of the rotated C 0 O 4 tetrahedra. The total 
energy change in this case is mainly controlled by relatively large anisotropy parameter Ki. 

In practice, such a situation can be realized by applying the magnetic field along the 2 : 
axis, which cants the spins out of the xy plane (Fig. (HI). The angle 0 between the spins and 
the xy plane can be estimated as cos^ = —y^gHSKSJ + 2i^i). There is a clear similarity 
with the atomic picture, depicted in Fig. [71 the magnetic field slightly decreases \Pz\ and 
induces two perpendicular components of the the polarization, which satisfy the condition 
Py = —Px‘ According to the atomic calculations (Fig. (T]), the components Py and Px are 
mainly induced at the Co sites 1 and 2, respectively, that is closely related to the geometry 
of the rotated C 0 O 4 tetrahedra. 
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FIG. 8 . (Color online) Electric polarization in the external magnetic field parallel to the 2 axis, as 
obtained from the solution of the low-energy electron model in the Hartree-Fock approximation. 
The starting configuration, where antiferromagnetic spins were parallel to the [110] axis, is shown 
in the inset. 

IV. SUMMARY AND CONCLUSIONS 

Using effective low-energy model, derived from the first-principles electronic structure 
calculations, we have investigated the multiferroic properties of Ba 2 CoGe 207 . There are two 
important aspects, which make this material especially interesting in the field of multiferroic 
applications: (i) Ba 2 CoGe 207 crystallizes in the noncentrosymmetric but nonpolar structure. 
The key symmetry operation, which controls the multiferroic properties of this material is the 
rotoinversion. From the view point of magnetic symmetry breaking, there is a fundamental 
difference between inversion and rotoinversion. The magnetic pattern, which can break the 
inversion symmetry, should be rather nontrivial. Typically, such magnetic order arises from 
the competition of many magnetic interactions and crucially depends on a delicate balance 
between these interactions. On the contrary, the rotoinversion symmetry can be broken 
very easily, simply by caning the magnetic moments out of the rotoinversion axis, (ii) The 
magnetic structure of Ba 2 CoGe 207 is very soft and can be easily deformed by the external 
magnetic field. This property is related to the specific geometry of the Ba 2 CoGe 207 lattice, 
where the C 0 O 4 tetrahedra are interconnected with each other via the Ge 04 octahedra. 
Thus, the magnetic Co atoms are separated by the long Co-O-Ge-O-Co paths, resulting 
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in the relatively weak exchange coupling J. Another important ingredient is the weak in¬ 
plane anisotropy, which is inherent to magnetic compounds with the spin 3/2.— We propose 
that this anisotropy has an intraatomic nature and related to the anisotropy of the orbital 
magnetization. 

Nevertheless, it seems that the magnetic softness of Ba 2 CoGe 207 has also one negative 
aspect: the Neel temperature Tn = 6.7 K, below which the multiferroic behavior has been 
observed, is relatively small.- This imposes a serious constraint on the practical realization of 
the considered effects: since Tn is controlled by the same exchange coupling J, any attempts 
to increase Tn will make it more difficult to deform the magnetic structure and manipulate 
the properties of Ba 2 CoGe 207 by the magnetic field. 

On the theoretical side, we have shown that the electric polarization of Ba 2 CoGe 207 can 
be presented as the sum of electric dipoles, which are induced at each Co site by the local 
exchange field. Each such dipole is given by the expectation value (p) = Tv[pD\ of the 
dipole matrix p and the site-diagonal density matrix T>. This is rather general property of 
Ba 2 CoGe 207 , which was derived starting from the Berry-phase theory of electric polarization 
in the Wannier basis. The local character of polarization in the case of Ba 2 CoGe 207 is 
directly related to the rotoinversion symmetry. 

The magnetic state dependence of the electric polarization is fully described by the site- 
diagonal density matrix P. Any rotation of the local magnetization from the rotoinversion 
axis lowers the symmetry of the density matrix T> at the Co site and induces the local electric 
dipole due to the transfer of the weight of the Wannier functions to the neighboring oxygen 
sites. This transfer is possible due to the p-d hybridization. However, the spin dependence of 
the hybridization itself does not play an important role. The direction of the electric dipole 
depends, via P, on that of the local magnetization. The total polarization of the crystal is 
the macroscopic average over the microscopic electric dipoles. This is the basic microscopic 
picture underlying the multiferroic behavior of Ba 2 CoGe 207 . 

Acknowledgements. This work is partly supported by the grant of Russian Science Foun¬ 
dation (project No. 14-12-00306). 
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Appendix: Polarization of oxygen band and interatomic exchange interactions 


In this appendix, we evaluate the change of the magnetic energy, caused by the polariza¬ 
tion of the O 2p band. 

After the solution of the low-energy model, consisting of the Co 3d bands, we expand the 
basis Wannier functions of the model in the original LMTO basis, Eq. ([TD]), and construct 
the spin magnetization density m(r) = n|(r)—associated with the Co 3d band. In this 
appendix, and (n-j- and ^’|) denote the electron densities (potentials) for the majority 
(t) and minority (4-) spin states. m(r) has major contributions at the Co sites as well as 
some hybridization-induced contribution at the oxygen and other atomic sites. Following the 
philosophy of the low-energy model,— the interaction of m(r) with the rest of the electronic 
states can be described in the frameworks of the local-spin-density approximation (LSDA). 
Therefore, our strategy is to evaluate, within LSDA, the exchange-correlation (xc) field 
h[r) = t;|(r)—n'|-(r), which is induced by m(r) and polarizes the O 2p band, and find the 
self-consistent change of m(r) and b(r), caused by the polarization of the O 2p band. For 
these purposes, it is convenient to use the SCLR theory.— For simplicity, let us consider 
the discrete lattice model and assume that all weights of m(r) are concentrated in the 
lattice points: m(r) = ~ where rriy is the local magnetic moment at the 

site V. Furthermore, we recall that LSDA is conceptually close to the Stoner model.— 
Then, the magnetic part of the xc energy can be approximated as E^c = ~ \ In 

practical calculations, the parameters {/i;} can be obtained using the values of intraatomic 
spin splitting and local magnetic moments. In LMTO, the intraatomic exchange splitting 
can be conveniently expressed via C-parameters of the centers of gravity of the Co 3(i states.— 

Then, by introducing the vector notations b = [by], and the tensors X = [Iy6yy>] and 
'JZ = [JZyy>], the self-consistent field can be found as 


b = 


l+X7^ 


1 -1 


b°, 


where b^ = Xm is the xc field induced by the Co 3d band, and the response tensor is 
obtained in the first order of the perturbation theory for the wavefunctions, starting from 
the nonmagnetic band structure in LDA: 


occ unocc BZ 


= y] y] y] y] 


(^va\*^va //^v'b\*^v'b 
^n'kV^n'k/ ^nk 


ab n 


^nk ^n'k 


+ C. C. 


(A.l) 
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In these notations, {C^^} coefficients of the expansion of the LDA wave functions 

over the LMTO basis, {^nk} are the LDA eigenvalues, and k runs over the first Brillouin 
zone (BZ). Moreover, similar to constrained random-phase approximation for the screened 
Coulomb interactions (Ref. Il8l ). we have to exclude from the summation in Eq. (lA.ip the 
contributions, where both indexes n and n' belong to the Co 3d band. In the present 
perturbation theory, such terms describe the change of the magnetization, which is caused 
by the LSDA potential in the Co 3d band. However, these effects are already taken into 
account in the low-energy model, where the LSDA part is replaced by a more rigorous 
unrestricted HF approximation with the screened Coulomb interactions. Therefore, such 
terms should be excluded at the level of SCLR calculations for the LSDA part. In practical 
calculations, n runs over the occupied O 2p bands and n' runs over the unoccupied Co 3d 
bands. 

Once b is known, the change of m, caused by the polarization of the oxygen band, can 
be found as 

dm = —iZb 


and corresponding change of the xc-field is Sb = XSm. Since O 2p band is fully occupied, the 
net change of the magnetic moment is identically equal to zero: Sm^ = 0, irrespectively 

on the type of the magnetic order. Nevertheless, the individual moments 6 m^ can be finite 
and contribute to the energy. The corresponding correction to the total energy consists of 
two parts: SE = SEco-o + ^Eq, where SEqo-o = —\bnXXrh is the interaction of with 
the “external” xc field, created by the Co 3d band, and SEq is the energy change caused 
by dm in the O 2p band. It also consists of two parts: SEq = SEgp + dEjc, where 5Esp is 
the single-particle energy, which can be found in the second order of Sb as SEsp = jSP^iZb^ 
SEdc = \SmFXSrh is the double-counting energy, and SmF is the row vector, corresponding 
to the column vector dm. Meanwhile, it is assumed that the magnetic energy of the Co 3d 
band itself is described by the low-energy model in the HF approximation. 

SE may have different values in the case of the FM and AFM alignment of spins in the 
xy plane. This difference additionally contribute to interatomic exchange interactions in the 
plane. 

In the P42im structure of Ba 2 CoGe 207 , there are three types of oxygen atoms: 01, 02, 
and 03, which are located in the Wyckoff positions 2c, 4e, and 8/, respectively.- The ob¬ 
tained parameters {A;} are 0.98, 1.92, 1.10, and 2.01 eV for Co, 01, 02, and 03, respectively. 
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The magnetic moments are listed in Table U and the energies are in Table ITII As expected, 


TABLE I. Local magnetic moments rriy, derived from the low-energy model for the Co 3d band, 
and the moments Sruy, caused by the polarization of the 0 2p band in the ferromagnetic (F) and 
C-type antiferromagnetic (C) state of Ba 2 CoGe 207 . All magnetic moments are in /iB per site and 
the number of such sites in the unit cell is given in the parentheses. 




F 



C 


rtiv 


Smy 

rriy 


Smy 

Co (xl) 

2.251 


0.304 

2.245 


0.234 

01 (xl) 

0.005 


-0.004 

0 


0 

02 (x2) 

0.004 


-0.004 

0 


0 

03 (x4) 

0.184 


-0.073 

0.180 


-0.057 


TABLE 11. Magnetic contributions to the energy of interaction between Co 3d and O 2p bands 
((5Eco-o)’ fhe single particle energy in the O 2p band {SE^p), the total energy of the 2p band 
(SEq), and the total energy {SE = SEqo-o + ^T'o) obtained for the ferromagnetic (F) and 
C-type antiferromagnetic (C) states. All energies are in meV per one formula unit. 



F 

C 

O 

o 

O 

-282.98 

-217.67 

6Esp 

-6.46 

-3.24 

6Eq 

27.07 

16.91 

6E 

-255.91 

-200.76 


the moments rriy are distributed mainly between central Co site and its neighboring sites 03. 
In the FM state, the total moment is mco + ttiqi + 2mo2 + 4mo3 = 3 /iB, which is totally 
consistent with the value obtained in the Wannier basis. The moments rriy and Srriy are 
parallel at the Co sites and antiparallel at the oxygen sites. This tendency is consistent with 
results of the first-principles calculations and can be explained by the hybridization between 
Co 3d and O 2p states.— Therefore, the negative sign of 6 Eco-o is due to the contributions 
of the Co sites, which are partly compensated by the positive contributions of the oxygen 
sites. The absolute value of SEco-o is larger for the FM state, mainly because mco and 
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Stuqq are larger than those in the C-type AFM state. Thus, the Co-0 interaction addition¬ 
ally stabilizes the FM order. Then, since the O 2p band is fully occupied, the contribution 
6 Esp is relatively small. The total contribution of the O 2p band to the magnetic energy is 
positive. This is because the fully occupied O 2p band itself does not favor the magnetism 
and any magnetic polarization of this band increases the energy. This also explains why 
5Eq is slightly smaller in the C-type AFM state: the magnetic moments 5m^ are smaller 
and, therefore, the magnetic perturbation of the O 2 p band is also smaller. Anyway, this 
effect is considerably weaker in comparison with the change of 5Eqq_q. 


In total, the magnetic polarization of the oxygen band favors the FM alignment of spins 
in the xy plane. By mapping the energy change 5E onto the spin model, which includes only 
nearest-neighbor interactions in the xy plane, the change of the exchange coupling, caused 
by the polarization of the oxygen band, can be estimated as SJ = |((5_E[C] — (5-E[F]) ^ 
14 meV. Thus, 6 J will indeed compensate the AFM exchange coupling, obtained in the 
low-energy model for isolated Co 3d bands. However, despite the correct tendency, the 
obtained change SJ is too large (and would lead to the FM alignment in the xy plane). 
There may be several reasons for it: (i) The SCLR theory may be to crude for treating the 
magnetic polarization of the oxygen band (in fact, the perturbation, which is introduced 
hy in the O 2p band is not small); (ii) Some quantitative estimates may change by 
considering the correct crystal structure below Tn (which is not available yet). Particularly, 
the original P42im symmetry can be lowered by the exchange striction effects; (iii) Tp some 
extent, the correlation interactions in the Co 3d band, beyond the HF approximation, will 
additionally stabilize the C-type AFM order. The corresponding contribution to the total 


energy difference between FM anc 


C-type AFM states can be estimated using the second- 


order perturbation theory (Ref. IlSl) as 2 meV per one formula unit. This change alone does 


not significantly change the FM contribution SJ, caused by the polarization of the O 2p 
band. However, it can collaborate with other effects, such as the exchange striction. 
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